Accurate electronic band gap of pure and functionalized graphane from GW 

calculations 
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Using the GW approximation, we study the electronic structure of the recently synthesized hy- 
drogenated graphene, named graphane. For both conformations, the minimum band gap is found 
to be direct at the F point, and it has a value of 5.4 eV in the stable chair conformation, where 
H atoms attach C atoms alternatively on opposite sides of the two dimensional carbon network. 
In the meta-stable boat conformation the energy gap is 4.9 eV. Fhen, using a supercell approach, 
the electronic structure of graphane was modified by introducing either an hydroxyl group or an 
H vacancy. In this last case, an impurity state appears at about 2 eV above the valence band 
maximum. 

PACS numbers: 81.05.Uw, 71.15.Mb, 71.10.-w 
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After the discovery of graphene [l| and of its extraordi- 
nary electronic properties 0,0, 0| the chemical functional- 
ization of graphene has become the focus of special inter- 
est in contemporary materials science. Being first a truly 
two-dimensional crystal and demonstrating high elec- 
tron mobility [l|], graphene is ideal for modern electron- 
ics which is essentially two-dimensionalQ- At the same 
time, the peculiar gapless "ultrarelativistic" energy spec- 
trum of graphene [3J, l4J makes the creation of a "carbon 
transistor" based onp—n junctions highly nontrivial, due 
to the anomalous transparency ("Klein tunneling" Q ) . 
Nanoscale graphene single-electron transistor has been 
demonstrated already Q but it has a relatively restricted 
domain of potential applications. Therefore, to trans- 
form graphene into a semiconductor with a conventional 
electron spectrum keeping its two-dimensionality is a real 
challenge in the field, and chemical functionalization is 
considered as one of the most promising ways to solve 
the problem. 

A new derivative of graphene was recently 
synthesized^, where hydrogenation turned graphene 
into what now is called graphane. Previously this 
compound was studied theoretically, and it is actually 
an excellent example of a new and interesting material 
which was predicted from first principles theory 0, 0| 
before it was synthesised experimentally. In the cal- 
culations of the electronic structure it was found that 
graphane is a semiconductor with a rather wide energy 
gap[E E|. In the experiment 0] the adsorption of 
hydrogen on graphene was indeed observed to result in 
that a gap opened up in the electron states. Hence the 
adsorption of hydrogen turned the highly conductive 
graphene into insulating graphane, in accordance with 
the theoretical predictions, [g, However, the exact 



value of the band gap is still experimentally unknown. 
Previous theoretical calculations giving values from 3.5 
eV[8] to 3.8 eV[10] use a standard density functional 
with generalized gradient approximation (GGA). It is 
well known that, in general, this approach is not reliable 
to calculate energy gaps in semiconductors often giving 
inaccurate results (the so called "gap problem" ) . 
At the same time, knowledge about the size of this 
gap is crucial in order to assess the possibility of using 
graphane in electronics applications and for use in 
devices. 

The reason for a band gap opening up when hydrogen 
is adsorbed on graphene is that sp 2 bonded C atoms be- 
come sp 3 bonded atoms, where three of the four covalent 
sp 3 bonds are saturated by C atoms and the fourth cova- 
lent bond is saturated by an H atom. This is for instance 
illustrated in the theoretical part of Ref|7| and is also 
known from the adsorption of single hydrogen atoms on 
graphene [l2j • In accordance to expectations this change 
in chemical binding is also reflected in the electronic 
structure, where sp 3 bonded C atoms do not display a 
so called it* peak in the x-ray absorption spectrum [F2|. 
Hence the basic reason for the drastic change in conduc- 
tivity is in accordance to expectations for sp 3 bonded car- 
bon (i.e. diamond), with an insulating behaviour. There 
is of course no reason to expect that the size of the band 
gap of graphane should be similar to that of diamond, 
even though there are similarities in the nature of the 
chemical binding. Hence measurements or accurate the- 
oretical calculations are needed, and in this article we 
present a calculation based on the GW approximation. 

Due to the well-known deficiencies of DFT[3| to treat 
excited states, (for example the band gaps in the LDA 
or GGA approximations are much smaller than the ex- 
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perimental values), we have used the GW approximation 
(GWA) of Hedin FJ, |l5[ to study the electronic struc- 
ture of graphanc. In this formalism, the Kohn-Sham 
equations[13[ are replaced by the quasiparticle equation: 
(T + V ext + Wkn(r) + J dV£(r, r', £ n (k))Vw(r') 

= £ n (k)V>kn(r) 

where T is the free-electron kinetic energy operator, V ex t 
the external potential due to the ion cores, Vh the Hartree 
potential, £ the self-energy operator, and E n (k) and 
V'kra(r) are respectively the quasiparticle energy and wave 
function. 

An adequate approximation for the self-energy opera- 
tor is to write it as the product of the Green's function 
and the screened Coulomb interaction W, which yields 
to the so-called GW approximation. 

£(r,r',a;) = J dw'G(r, r', u + u')e i5w 'W{r, r', to') 

Then, since the difference £ — V xc between the self-energy 
and the Kohn-Sham exchange and correlation potential 
is small, a perturbation theory is used to write the QP 
Hamiltonian: 

H qp = H KS + V xc ) 

Finally, the QP energies arc obtained by expanding the 
real part of self-energy to first order around e® FT (k) and 
assuming that the the QP wave function and Kohn- 
Sham wave function ^kn are identical: 
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(k) + Z„ k x 



[(*£f T |ReS(r,r', £n (k))|^) - <*£f T \V X D / T (r)\^f T 
where the QP renormalization factor Z„k is given by 



7 DFT 



= [1 - «f T ||;ReI](r,r', e „(k))|^)]-i. 



Therefore, it appears that the DFT eigenvalues e„(k) are 
corrected by the GW approximation, giving a practical 
scheme to compute reliable excited state properties [3, 



17, 1 



Here we have used the code VASP[19j (Vienna Ab- 
initio simulation package), implementing the projector 
augmented waves (PAW) method [2(| to compute the 
ground state and excited state properties of graphane. 
First, we obtained a reliable geometry of the structure by 
optimizing it for the two possible conformations (see be- 
low) , using the generalized gradient approximation 21 1 . 
Contrary to Sofo et alJ8|, we have chosen to simulate 
completely isolated sheets of graphane, although the re- 
sulting bandstructure is very close to that other Q. For 
all the calculations, we have used the default cut-off for 
the wavefunction. During the optimization of the struc- 
tures, a k-point grid 22] of 10 x 10 x 1 was used. Then, 



using the relaxed geometries, a final run was performed 
with a 16 x 16 x 1 grid. Finally, these ground state calcula- 
tions were used to compute the quasiparticle bandstruc- 
ture with VASP, following the method described in Ref. 
23l . Two hundred bands were used for the summation 
over bands in the calculation of the polarisability and the 
self-energy, and a cut-off of 150 eV was used for the size 
of the polarisability matrices. For details about the im- 
plementation of the GW app roximation within the PAW 



formalism, see also Refs. |2J, [25j. Notice that the GW ap 



proximation was applied successfully to study the prop- 
erties of graphene nanoribbons 2(| 27 1, of graphene 28 1, 
and of carbon nanotubes[29j]. 

Graphane is made of a hexagonal network of carbon 
atoms, with the orbitals being of the sp 3 kind and bonded 
to hydrogen atoms, with a carbon/hydrogen ratio of one. 
It exists in two conformersQ, the chair- like conformer 
and the boat like conformer. In the chair conformer, all 
C-C bonds are equivalent and the hydrogen atoms are 
bonded alternatively on each side of the carbon plane 
(see Fig.l of RefJH). The boat conformer has two types 
of C-C bonds, depending on whether the carbons atoms 
are bonded to hydrogen atoms which are on the same 
side of the carbon plane or not. Our calculations indi- 
cate that the chair-like conformer is favored by about 
0.051 eV/atom, in very good agreement with the value 
of Sofo et al.[l|, who found 0.055 eV/atom. Also, for the 
chair-like conformer the bond lengths are 1.53 A (for the 
C-C bonds) and 1.11 A (for the C-H bonds), while for 
the boat-like conformer, the C-C bonds are 1.57 A and 
1.54 A long, and the C-H bonds are 1.10 A long. These 
values are essentialy identical to the one of Ref. In 
FigQ] we show the partial density of states (pDOS) for 
]the stable chair conformer, and it is clear that the occu- 
pied H states are located over an energy interval starting 
from the valence band maximum and extending some 10 
eV below. As regards to the unoccupied states, the H 
and C states are also extended and cover a wide energy 
range. 



Being confident about the reliability of our ground 
state calculations, we now turn to the calculation of ex- 
cited state properties. The GGA (full lines) and GW 
(dots) bandstructures for both conformers are presented 
in Fig [2] and [3] Firstly, our calculated GGA bandstruc- 
ture for the chair conformation agrees well with the one 
of Sofo et al (notice that the bandstructure of the boat 
conformation was not reported in this paper). The mini- 
mum band gap is direct at the high symmetry point T and 
has a value of 3.5 eV. However, this value is dramatically 
changed when using the more reliable GW approxima- 
tion, in this case the value is 5.4 eV. The transitions at 
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FIG. 1: (Color online) The total (in black) and partial (in 
red for carbon and in blue for hydrogen) density of states of 
graphane in the chair conformation. The top of the valence 
bands is chosen as the zero energy. 




FIG. 2: (Color online) The GGA bandstructure (full lines) 
and the GW bandstructure (red dots) of graphane in the chair 
conformation. The top of the valence bands is chosen as the 
zero energy. 



the high symmetry points M and K are also significantly 
increased, from 10.8 eV and 12.2 eV with the GGA ap- 
proximation, to 13.7 eV and 15.9 eV with the GW ap- 
proximation (see Table HJ. For the boat conformation 
(Fig. [3]), the band gap is also found to be direct at the 
r point, with a value given by the GW approximation 
of 4.9 eV, whereas the GGA gives only 3.3 eV. The val- 
ues of the transitions at the X, S, and Y high-symmetry 
points are also notably corrected, see Table |U Our cal- 
culations show that graphane, in both conformers, can 
be qualified as a large band gap insulator, and the elec- 
tronic properties are expected to be quite different from 
the ones of graphene. This is qualitatively in agreement 
with observations. 

Since graphane is a wide band-gap material it is rel- 
evant to ask if defect states can be introduced easily, 
and manipulated for electronics applications. We have 
here considered two examples of defects in graphane, in 
the chair conformation: replacement of a H atom by an 




FIG. 3: (Color online) The GGA bandstructure (full lines) 
and the GW bandstructure (red dots) of graphane in the boat 
conformation. The top of the valence bands is chosen as the 
zero energy. 



TABLE I: Values in eV of the transition energies of graphane 
at some high-symmetry points of the Brillouin zone for both 
conformers (chair or boat). The minimum band gap occurs at 
the r point. The last two lines refer to calculations performed 
with a 2 x 2 supercell, in which either an hydroxyl group (OH) 
or a H vacancy was introduced. 



Conformation 


Transition 


GGA value (eV) 


GW value (eV) 


Chair 


► E c 


3.5 


5.4 




M„ -> M c 


10.8 


13.7 




K„ -> K c 


12.2 


15.9 


Boat 


T-u > Tc 


3.3 


5.1 




X v — > X c 


7.0 


9.0 




Su ■ + Sc 


10.7 


13.9 




Y — ► Y 

J- v r J- c 


9.4 


12.6 


Chair+ OH 




3.3 


5.0 


Chair+ H vacancy 




3.7 


5.4 



hydroxyl group (OH), and H vacancies. We have cal- 
culated the electronic structure of such systems using a 
3x3 supercell, taking fully into account the relaxation 
of the structure due to the hydroxyl group or the H va- 
cancy. For both cases, a 2 x 2 x 1 kpoints mesh was 
used to obtain reliable geometry, and a final calculation 
with 8x8x1 kpoints mesh was conducted to obtain well 
converged density of states. 

In Fig. [4j we present our density of states for 
graphane+OH. The partial density of states of oxygen 
is presented for 's' and 'p'states. It appears that the 
electronic structure of graphane is not significantly mod- 
ified by the replacement of a H atom by the hydroxyl 
group. While the O-s states are very widely distributed 
in energy, the O-p states appear mainly as a peak close 
to the maximum of the valence band, between —1 and 
—2 eV. Among the extra electrons brought by the oxy- 
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FIG. 4: (Color online) The total (in black) and partial (in 
red for O-s states and in blue for O-p states) density of states 
of graphane+OH. The top of the valence bands is chosen as 
the zero energy. 
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FIG. 5: The total density of states of graphane+ a H vacancy. 
The Fermi level is put at zero energy. 



gen atom, one electron is involved in a bond with the 
neighbouring C atom, another electron is participating 
in the bond with H, and the two remaining electrons are 
forming a doublet on the oxygen atom. This is confirmed 
by the fact that the geometry of the network of carbon 
atoms, driven by sp 3 hybridization, is almost unaffected 
by the addition of the hydroxyl group. The band gap 
(from GGA) is slightly decreased, from 3.5 eV for pure 
graphane, to 3.3 eV here. Here as well, the band gap 
is underestimated by GGA, so we have applied the GW 
approximation to get a meaningful value. However, since 
the numerical cost of the GW approximation is quite 
high, we had to reduce the size of the supercell to 2 x 2. 
In this case, the GW band gap is about 5.0 eV (see Tab|T| 
which is in line with the band gap of pure graphane (5.4 
eV). 

Next we introduced a hydrogen vacancy in graphane. 
Then, not every C atom of the graphane layer of the 
graphane layer is saturated with a H atom. From inspec- 
tion of the relaxed geometry, it appears that the bonds 
between the carbon neighbouring the H vacancy and the 
carbon network have gained a significant part of sp 2 hy- 
bridization, like in graphene, with bond lengths here of 
1.49 A (1.42 A in graphene). The resulting total density 
of states is shown in Fig|5] An impurity state develops 



close to 2 eV above the valence band maximum, at the 
Fermi level, changing significantly the overall picture of 
the DOS. Also, using the GW approximation (for a 2 x 2 
supercell), the difference in energy between the valence 
band maximum and the conduction band mininum is in- 
creased from 3.7 eV to 5.4 eV, but the impurity state at 
the Fermi level is not affected by the GW correction. 

To summarize, in this article, results concerning the 
electronic structure of graphane have been presented us- 
ing the GW approximation. It was found that for both 
conformers, the minimim band gap is direct and located 



at the r point. For the stable chair conformation, the 
most stable one, it has a value of 5.4 eV, while for the 
boat conformation, it has a value of 4.9 eV. We also found 
that defects in graphane, in the form of H vacancies, pro- 
vide an impurity state some 2 eV above the valence band 
maximum, a finding which may be important when try- 
ing to use this material for electronics applications. 
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